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Abstract Torrefaction is a thermal treatment step in a 
temperature range of 210-240 °C, which aims to improve 
the dimensional stability and durability of wood. The mass 
loss kinetics for torrefaction of wood samples was studied 
using equipment specially conceived to measure mass 
losses during thermal treatment. Laboratory experiments 
were performed under nitrogen for heating rates of 0.1, 
0.25, 1, and 2 °C min -1 . A mathematical model for the 
kinetics of the thermodegradation process was used and 
validated. Measurements of temperature distribution and 
anhydrous mass loss were performed on dry sample of 
poplar wood during torrefaction in an inert atmosphere for 
different temperatures. The mathematical formulation 
describing the simultaneous heat and mass transfers 
requires coupled nonlinear partial differential equations. 
These unsteady-state mathematical model equations were 
solved numerically by the commercial package FEMLAB 
for the temperature under different treatment conditions. A 
detailed discussion of the computational model and the 
solution algorithm is given below. Once the validity of 
different assumptions of the model had been analyzed, the 
experimental results were compared with those calculated 
by the model. Acceptable agreement was achieved. 
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Introduction 

Nowadays, the use of wood as building material is pro¬ 
moted. The heat treatment of the wood by mild pyrolysis is 
used to improve some of the characteristics of the final 
wood product, such as, its durability and dimensional sta¬ 
bility. This process substitutes the use of external chemical 
products, which enables the treated wood to remain as an 
environmentally friendly product. Its an effective method 
to improve biological durability of wood [1-3]. The heat 
treatment of the wood induces a chemical modification of 
the main wood constituents. The lignin polymer structure is 
modified [4-8], the ratio between the amorphous and the 
crystalline cellulose is also changed [9-11], hemicelluloses 
are strongly decomposed [10, 12, 13], and some precursor 
products of the charcoal appear [14]. These chemical 
modifications confer with the new wood properties, such 
as, the improved dimensional stability [15-17], the 
improved fungal resistance [18-20], some color changes 
[21-25], hydrophobic increase [19, 22, 25-29], and 
improved mechanical properties [30-32]. 

The concept of thermal treatment to stabilize the wood 
structure has lead to the development of several treatment 
processes in some European countries [1]. In the 70s, the 
research developed by the Ecole des Mines de Saint-Eti- 
enne (France) has worked out for the first time, an indus¬ 
trial process named «Retified wood®». More recently, in 
the 90s, the forest industry in Finland has developed a new 
process; Thermo wood®. Some others technologies have 
been developed, such as, the Netherlands pro¬ 
cess «PLATO®» and the Menz Holz’s process OHT (Oil 
Heat Treatment) in Germany. The wood torrefaction is a 
heat treatment of ligno-cellulosic material carried out at 
temperatures up to 300 °C. Torrefaction can be carried out 
under different gaseous atmospheres. Gaseous combustion 
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products released during the thermal treatment are 
employed most of the time. Chemical reactions involved 
during torrefaction, as well as final properties of the 
material depend strongly on the treatment temperature and 
the duration of the process. During torrefaction, slow 
pyrolysis prevails: wood is thermally decomposed at a slow 
rate [33]. The anhydrous mass loss (AML) resulting from 
the wood decomposition is representative of the physical- 
chemical transformations of wood. The AML matches to 
the advancement of the torrefaction process. 

According to the treatment intensity, the conferred 
properties of the wood vary strongly. For high values of the 
treatment temperature and time, the fungal resistance 
becomes total, the coloring is darkened, and the mechanical 
resistance is strongly distorted. Recent studies have shown 
the direct relation between the anhydrous mass loss (AML) 
due to thermodegradation and the conferred wood proper¬ 
ties [20-34]. Hence, controlling the heat treatment process 
means to control precisely the AML of the material. It 
requires a coupled model, taking into account, heat transfer 
and mass transfer due to thermodegradation reactions 
inside wood sample. 

The mechanism of wood thermodegradation shows the 
presence of several phases when the temperature increases. 
A phase of elimination of some volatile compounds takes 
place at temperatures lower than 190 °C. For temperature 
range of 190-280 °C, hemicelluloses are converted 
essentially into gases and acetic acid. This step corresponds 
to wood roasting. 

At the temperature range about 250-300 °C, the lignin 
and cellulose decompose to give three products: gas 
(vaporization of volatile extractives and of bound water 
absorbed on the wood fibers, formation of volatile by¬ 
products like furfural or hydroxymethylfurfural, the depo¬ 
lymerization of lignin network through cleavage linkage or 
dehydration reactions lead to the formation of numerous 
low molecular weight phenolic compounds), tar, and char. 
This decomposition leads to an anhydrous mass loss, rep¬ 
resentative for wood chemical transformation. Properties of 
heat-treated wood depend on the mass losses [35]. Hence, 
controlling the quality of the heat-treated wood means to 
control precisely the mass loss during the treatment. 
Therefore, the knowledge of the kinetic schemes of the 
wood thermodegradation is required. 

The aim of this study is to develop a tool allowing the 
prediction of the temperature profiles in the wood and the 
associated AML during the treatment. A two-stage semi- 
global kinetic model is used here to describe the wood 
thermodegradation. Numerous experiences of heat treat¬ 
ment at 240 °C under nitrogen were realized with poplar 
wood. The curves of the instantaneous mass loss and 
temperature were recorded. The experimental data were 
used to determine the model parameters. Finally, a 


comparison between the experimental data and the 
numerical results was realized. 


Materials and methods 

Wood sample 

Heat treatment is carried out on wood panels of large 
sections, dried in the steam room at 105 °C, until their 
mass stabilizes around 269-280 g. The panel’s dimensions 
are 25 x 11 x 2.5 cm 3 , respectively, in the longitudinal, 
tangent, and radial directions. The wood specie used in the 
present study is the Poplar (Populus nigra). 

Heat treatment by conduction and measure 
of the evolution of the dynamic sample mass 

Heat treatment has been carried out under nitrogen by 
conduction between two metallic heating plates placed on a 
precision balance (Fig. 1). The instantaneous mass and 
temperatures at the surface and in the middle of the sample 
are recorded. The heat treatment phase consists of three 
successive thermal stabilization areas (Fig. 2). The first one 
is kept at 105 °C to verify the mass stability (anhydrous 
mass). The second one (165 °C) is useful for the thermal 
homogenization of the panel. The last one corresponds to 
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Fig. 1 Schematic equipment used for heat treatment. 
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Fig. 2 Typical heat treatment schedule. 
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the thermodegradation operating temperature. Each test has 
been repeated three times to check the reproducibility. 


Mathematical formulation 


Modeling of the heat and mass transfers 


Wood is a complex material whose thermodegradation 
involves a series of chemical reactions with coupled heat and 
mass transfer. This section is devoted to the mathematical 
formulation of the mentioned transfers in the wood porous 
media introducion, as well as the heat effects due to the reac¬ 
tions of polymer decomposition, as the formation of volatiles 
and non-degradable solid products. The following simple 
assumptions are used for the mathematical formulation: 

- the wood sample has a rectangular form characterized 
by its length L, width / and thickness e , respectively in 
the longitudinal, tangent, and radial directions. The 
sample thickness is very weak compared with other 
dimension and the modeling can be represented in ID. 

- the convection transfer occurring because of the 
formation of volatiles inside the wood material during 
the heat treatment is neglected. 

- The dimensions of the sample are considered invari¬ 
able; the formation of fissures and cracks is neglected. 

- The wood physical and thermal properties, respec¬ 
tively, the specific heat Cp , and the wood thermal 
conductivity X are considered as a linear function of the 
temperature [36]. 


Cp = 1112 + 4.85(7’ -273) 

(1) 

A = 0.13 + 0.0003 (T- 273) 

(2) 


The macroscopic conservation equation governing the 
heat transfer phenomenon is given by: 


P Cp 


dT 

dt 



dp 


+H *t 


(3) 


where p = m(t)/V is the averaged anhydrous wood density, 
T is the local temperature, and H p is reaction enthalpy. The 
instantaneous sample mass m(t) is detailed in Eq. 11, the 
sample volume V is considered invariable. 

The process of the wood heat treatment involves several 
reactions of polymer degradation and synthesis giving 
numerous reaction products [35, 37]. Consequently, mathe¬ 
matical modeling describing the wood thermodegradation is 
a simplified approximation of the real phenomena. However, 


the knowledge of kinetic schemes describing the thermode¬ 
gradation process is required for the elaboration of predicting 
tools that allow the control of wood heat treatment. 

Numerous mechanisms have been published to provide 
a rational explanation of the wood thermodegradation and 
the associated heat effects. Three main approaches could 
be distinguished: the first one considers the wood as a 
single homogeneous element [38-41], the second one takes 
into account the decomposition of the essential wood 
compounds: cellulose, hemicelluloses, and lignin [42], and 
the third one considers the existence of three fractions 
(pseudo-compounds) whose decomposition obeys different 
laws [43]. It is usually assumed that, the reaction kinetics 
follows an Arrhenius law with two parameters: activation 
energy and kinetic constant. 

The mechanisms adopted in this study, is based on a 
two-stage, semi-global, multi-reaction kinetic model of 
wood mild pyrolysis, where the wood is subdivided into 
three pseudo-components A h A 2 , and A 3 , as it was previ¬ 
ously proposed [43]. Each pseudo-component is charac¬ 
terized by a specific kinetic law of decomposition and its 
mass fraction, respectively a 1? a 2 , and a 3 , such as, 
a 3 = 1 — + a 2 . At a temperature level lower than 

200 °C, the elimination of some volatile wood compounds 
occur. At the temperature range of 165-240 °C, the deg¬ 
radation of Ai prevails, giving a gaseous product G x . The 
thermodegradation of A 2 leads to the production of a solid 
fraction C 2 and a gas product G 2 . The proposed kinetic 
scheme is given by the Eqs. 4-5. 

Ai—+*iG! (4) 

A 2 —* k 2 y 2 C 2 + (1 — y 2 )G 2 (5) 


The mass fractions of the non-degradable solid y 2 
depends on the temperature. 

The assumption that, the kinetics of all the involved 
reactions is described by the first order laws, and the mass 
balance equations can be written, respectively as follows: 


dm Al 

dt 

dm A2 

dt 


~k\m Al 


—k 2 m Al 


dm C2 

dt 


ki y 2 m Al 


( 6 ) 

(?) 

( 8 ) 


where m Al , m Al and mc 2 are respectively the mass of the 
constituents A 1? A 2 , and C 2 ; k x and k 2 are the reaction rate 
constants obeying the Arrhenius law such as: 


Table 1 The kinetic parameters used in the numerical simulation for the poplar. 


Species 

a i 

0( 2 

Ea i/kJ mol 1 

Ea 2 t kJ mol 1 

Ws 

* 0 l/s 

2 

H p /kJ mol 1 

Poplar 

0.08 

0.12 

118.5 

119 

4.10 7 

6.10 8 

0.07 

123 
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Fig. 3 Time evolution of average weight loss, surface temperature, 
and center temperature: Poplar (1 °C min -1 , 7 max = 210 °C). 

k\ = k 0 i exp (9) 

k 2 =ko2 ( 10 ) 

where, k 0] and k 02 are the pre-exponential factors; Eci\ and 
Ea 2 are the activation energies. The instantaneous total 
mass of the sample is equal to: 


Fig. 4 Time evolution of average weight loss, surface temperature, 
and center temperature: Poplar (1 °C min -1 , T max = 220 °C). 

m(t) = m Al (t) + m Al ( t) + m A ,(t) + m C2 {t) (11) 

Determination of the model’s parameters 

The resolution of the Eqs. 3-10 requires the determination 
of 8 parameters: the pre-exponential factors k 0 1 , k 02 , and 
activation energies Ea u Ea 2 , the mass fraction fractions a 1? 
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Fig. 5 Time evolution of average weight loss, surface temperature, 
and center temperature: Poplar (1 °C min -1 , T max = 230 °C). 

a 2 , and y 2 \ the reaction enthalpy Hp. These parameters 
have been obtained by minimizing an error function 
between the experimental data and simulation results. The 
number of parameters to determine is large and the opti¬ 
mization procedure becomes difficult. Consequently, it is 
necessary to find a good method to obtain an efficient 
optimization. Because of the complexity of the model, the 





Fig. 6 Time evolution of average weight loss, surface temperature, 
and center temperature: Poplar (0.25 °C min -1 , T max = 230 °C). 

dependence between the parameters, and the existence of 
several local optima, the deterministic methods (conjugate 
gradient, Lagrange etc.) could not be used. A method based 
on a random selection of numbers and exploring the search 
space more efficiently is required. As example of such kind 
of methods, it can be cited to the Monte Carlo algorithm or 
genetic revolution [44-46], the hybrid algorithm, the sim¬ 
ulated annealing algorithm, etc. The genetic revolution has 
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been adopted in this study. The error function to minimize 
over both the thermal and mass measurements is given in 
the Eq. 12. 

1 |^exp(0 ^calculated (0 | 

error global =~2^ (r exp (/) + r calculated (/)) 

|AML exp (/) AML ca i cu i atec i(/) | 
(AML exp (/) -f AML ca i cu ] atec i (/)) 

( 12 ) 

where n x is the number of data points of the temperature, 
r exp (/) is the experimental temperature, T calculated (z) is the 
calculated temperature on the i -th grid point, n 2 is the 
number of data points of the mass loss AML exp (/') and 
AML calculated (/') are respectively the experimental and the 
calculated mass loss due to the thermodegradation on the j- 
th grid point. The anhydrous mass lost is defined as 
follows: 

AWL = (m 0 -m(t))/m 0 (13) 

m 0 is the anhydrous wood mass before heat treatment. 

The experiments are carried out in such a way, allowing 
the determination of the model parameters over two stages. 
At the first, with the isothermal area at 165 °C, the reaction 
given by the Eq. 5 can be neglected, because of its very 
slow rate. It has also been considered that the reaction 
enthalpy of the first reaction is negligible. The three 
parameters relative to the first reaction (Eq. 4) are then 
determined: oq, k ox , and Ea x . The second isothermal area at 
240 °C, corresponding to the heat treatment process, 
allows the determination of the five other parameters: a 2 , 
y 2 , & 02 , Ea 2 , and Hp. A random algorithm has been used, 
and the minimum and maximum values of each parameter 
have been initialised according to the values available in 
the literature [39, 42, 43, 47]. 

The governing equations are solved numerically with 
the commercial package, FEMLAB [48]. This is designed 
to simulate systems of coupled, non-linear, and time 
dependent partial differential equations (PDE) in one-, two- 
, or three- dimensions. The mesh convergence was verified 
with refined mesh sizes. Time step of 360s, 1866 nodes, 
and 8080 elements mesh size were considered to be 
appropriate. 

The kinetic parameters obtained using the method 
described are in Table 1. 


Results and discussion 

Effect of maximum temperature 

This section is devoted to the investigation of the effect of 
temperature of the heat treatment on the anhydrous mass 





Fig. 7 Time evolution of average weight loss, surface temperature, 
and center temperature: Poplar (0.5 °C min -1 , r max = 230 °C). 

loss, the parameter, chosen as the controlling parameter for 
the heat treatment process. Three experiments have been 
carried out using poplar wood during 20 h under heat 
treatment at temperatures 210, 220, and 230 °C, respec¬ 
tively. The heat treatment d rate has been fixed to 
1 °C min -1 . Experimental and numerical results are shown 
on Figs. 3, 4, and 5. It can be observed that the increase of 
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Fig. 8 Time evolution of average weight loss, surface temperature, 
and center temperature: Poplar (2 °C min -1 , T max = 230 °C). 

the treatment temperature enhances the AML. A good 
agreement is observed between the experience and the 
numerical simulation. 

Effect of heating rate 

In this paragraph, the study has been performed on the 
poplar wood, under treatment temperature of 230 °C and 


for different heating rate of 0.25, 0.5, 1, and 2 °C min 1 , 
respectively. Results are shown on Figs. 5, 6, 7, and 8. 


Conclusions 

The kinetics for the pyrolysis of wood was investigated. 
Based on the laboratory experiments performed with the 
described device, a mathematical model for the kinetics of 
the pyrolysis process was proposed and validated. This 
model incorporates the reaction pathway that was devel¬ 
oped by Grioui et al. [43], for the description of the 
transformation of wood subject to heat. This model allows 
the computation of internal profiles of temperature and 
average mass loss. The Femlab [48] software was used to 
solve the system of partial differential equations. This 
study includes a validation section, that shows the ability of 
the model to accurately predict the internal profiles of 
temperature, and mass loss in the experimental configura¬ 
tion. Comparison of the simulated temperature profiles and 
average mass loss with experimentally measured values 
showed reasonable agreement. 

Further studies may allow more accurate modeling of 
wood torrefaction mainly, by taking into account the 
thermodegradation of lignin and accurate thermo physical 
properties. This model can be used in the design of 
industrial pyrolysis installations. 
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